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INTRODUCTION 


External  beam  radiation  therapy  (EBRT)  is  a  frequently  used  tool  in  the  treatment  and 
management  of  prostate  cancer.  One  of  the  standard  treatments  for  prostate  cancer  has  been  three 
dimensional  conformal  radiation  therapy  (3DCRT) ,  in  which  the  radiation  fields  are  conformed  to 
the  shape  of  the  prostate  at  each  angle,  allowing  for  irradiation  of  the  prostate  while  limiting  the 
irradiated  volume  of  the  surrounding  normal  tissues  [1,2].  Recently,  intensity  modulated  radiation 
therapy  (IMRT)  has  become  an  increasingly  popular  form  of  radiation  therapy  for  many  forms  of 
cancer,  including  prostate  cancer.  IMRT,  a  more  complex  radiation  delivery  scheme  than  3DCRT, 
uses  beams  which  are  modulated  throughout  treatment  to  produce  radiation  fields  of  non-uniform 
intensity  [3-5].  These  non-uniform  fields  provide  three  dimensional  dose  distributions  which  more 
closely  conform  to  the  shape  and  volume  of  the  prostate  than  3DCRT;  with  greater  conformity 
comes  a  larger  reduction  of  dose  to  the  surrounding  normal  tissues,  limiting  radiation  induced 
complications  of  treatment  and  organ  failure  [6-10]. 

The  greater  complexity  of  IMRT  results  in  longer  treatment  times  for  a  single  treatment 
session  or  fraction  of  radiation.  The  rising  prevalence  of  IMRT  as  a  treatment  modality  has  given 
rise  to  a  number  of  studies  investigating  the  impact  of  temporal  effects  in  IMRT,  especially  as 
they  relate  to  the  radiobiology  of  the  underlying  tissue  [11-17].  There  is  thus  an  increasing 
recognition  of  the  interaction  between  increasing  single  fraction  radiation  therapy  treatment  times 
and  cell  repair.  The  capacity  for  cell  repair  directly  effects  cell  survival,  thus  the  understanding 
of  these  temporal  effects  may  have  clinically  relevant  implications  with  regard  to  tumor  control 
and  normal  tissue  sparing  [18-24]. 

The  temporal  effects  relating  to  radiation  therapy  (RT)  are  generally  analyzed  in  terms  of 
cell  survival  using  the  linear  quadratic  (LQ)  model  [25-28].  For  a  single  IMRT  fraction,  the  LQ 
model  is  given  by 

5  =  exp-(«D-ty0G(O£>")  (1) 

where  D  is  the  dose  per  fraction  (dose/fx),  a  and  |3  relate  to  the  tissue’s  ability  to  repair  radiation 
damage,  and  G(t)  is  the  Lea-Catcheside  dose  protraction  factor,  also  known  simply  as  the  dose 
protraction  factor.  G(t)  accounts  for  the  time  variables  involved  in  RT.  In  its  most  general  form, 
G(t)  is  given,  at  a  time  t,  by  [29] : 


^  t  w 

^0  0 

where  D  is  the  dose  per  fraction  (dose/fx),  I(k)  is  the  dose  rate  at  time  k,  //  =  ln2/rj/2  ^  where  T1/2 

is  a  tissue  specific  repair  half-time,  and  w  and  v  are  time  variables.  The  LQ  model  is  dependent 
on  two  time-related  factors  which  correlate  to  radiobiological  effects.  The  first  of  these  is  the 
“fraction  duration,”  the  total  time  measured  beginning  when  the  first  radiation  beam  is  “turned 
on”  and  ending  when  the  last  beam  is  “turned  off.”  The  limits  of  integration  of  the  outer  integral 
in  Eq.  2  can  be  viewed  as  the  beginning  and  end  of  the  time  span  defined  by  the  fraction  duration 
(Tf).  The  effect  of  fraction  duration  on  cell  kill  has  been  well  studied,  with  an  increasing  fraction 
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duration  yielding  a  larger  surviving  fraetion  [11,  13,  16,  17,  19-24],  The  other  time  variable 
which  may  affect  the  degree  of  cell  kill  is  the  functional  form  of  how  the  dose  is  applied  over 
time,  or  “temporal  pattern  of  applied  dose.”  This  is  due,  in  Eq.  2,  to  the  double  integration. 
There  has  been  some  previous  implication  that  the  temporal  pattern  of  applied  dose  during 
treatment  may  impact  the  degree  of  cell  kill  [18,  27,  30-33].  However,  we  performed  the  first  in- 
depth  investigation  of  the  effect  of  the  temporal  pattern  of  applied  dose,  especially  with  regards 
toIMRT  [34]. 

This  model-based  study  revealed  techniques  for  temporal  manipulation  of  the  dose 
arrangement  to  maximize  or  minimize  cell  kill.  Cell  kill  was  found  to  be  maximized  using  a 
“Triangle”  dose  pattern  in  which  the  highest  dose  fields  of  a  plan  are  clustered  in  the  middle  of  a 
treatment  fraction  and  the  lowest  doses  distributed  between  the  beginning  and  end.  A  “V- 
shaped”  pattern,  in  which  the  lowest  dose  fields  were  clustered  at  the  center  of  a  treatment 
fraction  and  the  highest  dose  fields  are  distributed  between  the  beginning  and  end  of  a  treatment 
fraction,  was  found  to  minimize  cell  kill.  The  study  also  investigated  the  impact  of  three 
variables  associated  with  a  single  fraction  of  IMRT  that  can  interact  with  variations  in  the 
temporal  pattern  of  applied  dose  to  affect  cell  kill:  irradiated  tissue  type  (in  terms  of  cx/|3), 
fraction  duration,  and  dose  per  fraction  (dose/fx).  The  modeling  study  found  that  the  effects  of 
the  temporal  pattern  of  applied  dose  are  greatest  in  treatments  which  irradiate  tissues  with  lower 
oc/p,  use  longer  fraction  durations  and/or  large  doses/fx. 

The  “optimized”  dose  patterns  described  by  the  modeling  study  could  potentially 
influence  prostate  cancer  treatment  outcome.  This  is  because  integration  of  temporal 
optimization  a  treatment  planning  algortithm  could  be  used  to  enhance  tumor  cell  kill,  improving 
prostate  tumor  control,  or  decrease  cell  kill  in  normal  tissues  such  as  bladder  and  rectum  to 
reduce  radiation  sequeulae.  However,  before  conclusions  of  the  modeling  study  can  be 
integrated  into  clinical  techniques,  they  must  first  be  validated  in  vitro.  Then  steps  need  to  be 
taken  to  integrate  temporal  optimization  into  a  treatment  environment  and  test  the  ability  to 
delivery  such  a  treatment.  This  begins  with  the  development  of  a  phantom  for  use  in  three 
dimensional  radiobiology  experiments  in  a  realistic  treatment  environment.  Then  a  scheme  to 
implement  temporal  optimization  into  an  IMRT-based  prostate  cancer  treatment  plan  needs  to  be 
developed  and  analyzed. 

Our  hypothesis  is  that  the  temporal  pattern  of  applied  dose  in  a  single  fraction  of  radiation 
can  be  optimized  to  maximize  or  minimize  cell  kill.  Furthermore,  irradiation  with  a  treatment 
plan  incorporating  such  techniques  will  provide  differential  cell  kill  across  various  structures  of 
interest,  such  as  prostate  tumor,  bladder,  and  rectal  tissues,  compared  with  irradiation  using  a 
treatment  plan  derived  from  standard-based  methods. 

BODY 


There  are  three  main  tasks  associated  with  this  project.  Within  each  there  are  one  or  more 
subtasks,  each  of  with  has  several  objectives.  Progress  made  on  the  following  tasks 
programmed  to  be  addressed  in  the  first  year  of  the  project  is  as  described  below: 

Task  I:  Validate  the  modeling  study  in  vitro 
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A  manuscript  based  on  the  data  and  conclusions  resulting  from  accomplishment  of  the  objectives 
and  subtasks  of  Task  I  has  been  submitted  for  publication,  although  not  as  of  yet  accepted.  The 
discussion  of  the  accomplishment  of  these  subtasks  and  objectives  below  contains  excerpts  from 
this  submission. 

LA:  Experimentally  determine  radiobiological  parameters  of  cell  lines  of  interest  in  vitro.  (1-5 
mos.) 

Objectives: 

1.  Determine  cell  lines  of  interest.  (1-2  mos.) 

2.  Perform  in  vitro  experiments  on  all  cell  lines  of  interest  to  measure  parameters, 

repeating  experiments  to  ensure  result  consistency  as  needed.(2-4  mos.) 

3.  Analyze  survival  data  to  determine  radiobiological  parameters.  (4-5  mos.) 


As  noted  above,  the  modeling  study  predicted  that  temporal  optimization  effects  were 
greatest  for  tissues  with  low  a/|3.  For  a  comparative  analysis,  the  “cell  lines  of  interest”  as 
defined  in  Subtask  LA  Objective  1  should  include  at  least  one  low  cx/|3  line  and  one  high  a/|3 
line.  It  should  also  be  noted  that  the  modeling  study  held  a  constant,  meaning  there  was  also  an 
implicit  dependence  on  |3.  As  the  temporal  factors  couple  to  the  |3  term  in  the  LQ  model  (Eq.  1), 
the  impact  of  the  temporal  pattern  of  dose  will  be  stronger  in  cases  with  higher  (3.  However,  in 
terms  of  overall  cell  survival,  a  cell  line  more  susceptible  to  repairable  relative  to  irreparable 
damage  (thus  low  cx/(3)  is  also  important.  Ultimately,  this  means  that  the  greatest  opportunity  to 
see  differences  is  if  the  low  cx/|3  cell  lines  of  interest  for  this  study  also  have  high  (3,  while  the 
high  (x/(3  cell  lines  used  have  low  |3.  As  we  plan  to  measure  cell  survival  following  radiation 
using  the  clonogenic  cell  survival  assay,  we  also  require  that  the  cell  lines  we  use  be 
clonogenically  useful.  That  is,  they  should  form  good,  countable  colonies  for  clonogenic 
evaluation. 

Following  investigation  of  the  on-site  departmental  tissue  culture  bank,  we  found  three 
candidates  which  satisfied  our  requirements:  WiDr,  a  colonic  adenocarcinoma  line,  PC-3,  a 
prostate  carcinoma  line,  and  SQ-20B,  a  head-and-neck  squamous  cell  carcinoma  line.  Previously 
acquired  clonogenic  survival  data  for  these  cell  lines  implied  that  that  PC-3  and  WiDr  were 
lower  cx/|3  and  higher  (3  cell  lines,  while  SQ-20B  had  a  higher  cx/(3  and  lower  |3. 

Subtask  LA  Objectives  2-3  involved  the  experimental  measurement  of  the  radiobiological 
parameters  of  each  cell  line  used  in  the  completion  of  Task  1.  Measurement  of  a  and  (3  from  Eq. 
1  intertwines  with  the  previous  objective,  in  that  the  measurement  of  these  parameters  verifies 
that  they  indeed  satisfy  the  requirements  outlined  above,  a  and  |3  were  measured  by  creating 
survival  curves  for  each  cell  line  using  the  clonogenic  survival  assay.  The  survival  curves  were 
created  by  irradiating  cells  over  a  range  of  doses  from  0-10  Gy  and  determining  the  percentage  of 
cells  surviving  at  each  of  the  dose  points  analyzed.  Acute  irradiations  were  used,  eliminating  the 
impact  of  the  dose  protraction  factor  in  Eq.  1  (i.e.  as  the  dose  rate  ,  I  in  Eq.  2,  gets  very  large, 
G(t)  — >1).  Eq.  1,  with  G(t)=l,  is  then  fit  to  the  resultant  percent  survival  versus  acute  dose 
curves  to  determine  a  and  |3. 
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Cells  were  maintained  in  appropriate  media  at  37°  C  and  7%  CO2.  Taking  into 
consideration  the  cell  line  and  experimental  treatment,  cells  were  plated  in  6-well  plates  at 
densities  predetermined  to  yield  60-100  surviving  colonies  per  well.  A  non-irradiated  control 
with  100  cells/well  was  also  plated  for  each  cell  line  and  experiment  performed.  Prior  to 
irradiation,  cells  were  incubated  for  12-24  h  in  a  final  volume  of  2  ml/well. 

Cells  were  then  irradiated  with  each  cell  line  being  irradiated  individually.  After 
radiation  application,  cells  were  incubated  at  37°C  and  7%  CO2  for  10  days,  at  which  time  a 
clonogenic  survival  assay  was  performed.  Plates  were  washed  with  IX  phosphate  buffer 
solution,  and  then  stained  with  2%  crystal  violet  in  50%  ethanol  for  one  hour.  Excess  stain  was 
washed  off  and  surviving  colonies  (containing  >50  cells)  were  counted  by  eye.  The  treatment 
(Triangle  or  V-shaped)  of  the  plates  was  blinded  to  the  observer.  A  single  observer  counted  all 
events  and  comparisons  were  done  across  four  independent  observers  to  validate  the  results.  The 
average  number  of  events  per  plate  exceeded  400.  Surviving  fraction  was  calculated  for  each 
well  after  factoring  in  the  plating  efficiency  (PE)  determined  from  the  control  plates.  Each 
experiment  (i.e.  a  specific  cell  line,  fraction  duration  and  total  dose  combination)  was  performed 
multiple  times. 

Average  parameters  measured  using  these  techniques  validated  these  original  assertions 
regarding  each  cell  line,  with  WiDr  having  cx/(3=3.39  Gy  ((3=0.046  Gy'^),  PC-3  having  cx/(3=5.56 
Gy  ((3=0.046  Gy^),  and  SQ-20B  having  a/(3=8.72  Gy  ((3=0.016  Gy^).  As  “lower  a/(3”  cell  lines 
(WiDr  and  PC-3)  also  had  higher  (3  values  than  SQ-20B,  we  will  use  the  low  vs.  high  cx/(3  or  high 
vs.  low  (3  designations  interchangeably  in  this  report. 

The  measurement  of  the  radiobiological  parameters  related  to  the  G(t)  function  proved  far 
more  difficult.  Initially,  this  might  not  seem  the  case,  with  only  one  tissue  specific  parameter,  p, 
in  Eq.  2.  While  p  is  typically  defined  as  ln(2)/Ti/2,  Eowler  et  al.  [11]  notes  that  T1/2  continuously 
changes  with  time  and  that  two  T1/2  values,  with  properly  weighted  corresponding  G(t)  functions, 
is  a  reasonable  approximation.  This  is  the  approach  taken  in  the  previously  discussed  modeling 
study.  Thus  to  get  an  accurate  model,  at  least  four  additional  parameters  (two  T1/2S  and  two 
weights)  would  need  to  be  measured  for  each  cell  line.  A  series  of  “split  dose”  experiments,  in 
which  cells  were  irradiated  with  two  equal  acute  doses  split  by  a  time,  T,  were  performed  on  each 
cell  line.  In  this  specialized  case,  G(t)  reduces  to: 

G(0  =  ^(l  +  e-«")+^(l  +  e-^^")  (3) 

where  Wx  are  the  weights  and  px  are  the  corresponding  p  values.  Eig.  1  shows  an  example  PC-3 
split  dose  experiments  over  a  range  of  T  values  from  0-40  min.  Recall  that  the  surviving  fraction 
can  be  related  to  G(t)  from  Eq.  1.  Using  the  measured  a  and  (3  values  for  a  given  cell  line,  we 
attempted  to  fit  Eq.  3  to  the  data  measured  for  each  cell  line,  but  given  the  errors  in  each  data 
point,  we  found  the  parameter  fitting  would  either  give  unrealistic  results  or  confidence  intervals 
so  wide  that  virtually  any  measured  data  would  fit  the  model.  One  option  to  resolve  this  would 
perhaps  be  to  acquire  date  with  a  finer  split  time  sampling  than  we  have  done.  However,  the 
inherent  errors  of  the  clonogenic  assay  means  that  much  finer  time  sampling  would  not  give 
useful  results:  for  example,  in  Eig.l,  calculations  with  the  Mann-Whitney  U  test  show  the 
differences  in  surviving  fraction  between  all  points  >10  minutes  split  time  do  not  show  statistical 
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significance  (p  >  0.05).  Other  methods  we  used  to  attempt  to  estimate  G(t)  from  experimental 
data  for  more  complex  temporal  dose  patterns  suffered  similarly  from  the  inability  to 
differentiate  temporal  data  points  using  the  clonogenic  assay. 

We  were  thus  unable  to  directly  measure  all  the  radiobiological  parameters  to  determine 
if  the  temporal  optimization  effects  analyzed  below  in  Subtask  II.B  matched  the  radiobiological 
models  of  Eqs.  1-2.  Our  alternative  solution  was  to  analyze  the  entire  parameter  space  of  the 
modeling  study  to  validate  its  conclusions.  In  other  words,  will  we  strive,  as  originally 
programmed,  to  show  that  there  is  a  statistically  significant  difference  in  cell  survival  between 
irradiation  with  a  Triangle  versus  V-shaped  dose  pattern  for  low  0c/(3  cell  lines  at  high  dose/fx 
and  long  Tf.  Then,  however,  we  will  analyze  the  comparative  effects  of  temporal  optimization 
using  a  higher  cx/p  cell  line,  low  dose/fx,  or  short  Tf. 


PCS  Split  Dose  Experiment 


Split  Time  (Min) 

Figure  1:  Results  from  an  example  split  dose  experiment  on  PC-3  cells.  Error  bars  represent  the  standard  error 
at  each  analyzed  split  time.  Total  dose  was  6  Gy  for  all  points. 


LB:  In  vitro  experimental  verification  of  previous  modeling  study.  (3-6  mos.) 

Objectives: 

1 .  Perform  in  vitro  experiments  on  cell  lines  of  interest,  repeating  experiments  to  ensure 

result  consistency  as  needed.  (3-5  mos  ) 

2.  Analyze  of  survival  data  and  assess  results  in  comparison  to  assertions  of  the  modeling 
study.  (5-6  mos.) 
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WiDr,  PC-3,  and  SQ-20B  cells  were  all  irradiated  and  survival  evaluated  in  Subtask  LB 
using  the  clonogenic  assay  as  described  above  in  Subtask  LA.  As  the  dose  patterns  in  Subtask 
LB  are  designed  to  emulate  clinically  used  IMRT  irradiation  schemes,  irradiation  was  performed 
using  a  clinical  linear  accelerator  (2100CD,  Varian  Medical  Systems,  Palo  Alto,  CA)  delivering  6 
MV  photons  at  600  MU/min.  Each  six-well  plate  was  irradiated  individually.  The  plate  was 
placed  in  a  phantom  composed  of  a  rectangular  block  of  solid  water  with  a  plate-sized  cavity  (8.6 
X  12.8  X  2.0  cm)  at  the  center.  Radiation  was  delivered  with  a  20  x  20  cm  field  at  90  cm  SSD 
and  10  cm  depth  (i.e.  to  the  bottom  of  the  plate,  where  the  cells  were).  Dosimetry  in  the  phantom 
was  assessed  with  film  and  monitor  units  were  calculated  accordingly  to  deliver  the  desired  dose. 

Once  in  the  phantom,  each  plate  was  irradiated  with  900  cGy  divided  into  six  fields, 
delivered  over  20  minutes  (Tf  =20  min).  The  six  fields  were  comprised  of  two  fields  each  of  325 
cGy,  100  cGy,  and  25  cGy  with  the  temporal  midpoints  (times  at  which  half  the  dose  of  an 
individual  field  has  been  delivered)  equally  distributed  over  Tf.  Two  plates  of  each  cell  line  were 
irradiated,  one  with  the  fields  arranged  in  a  Triangle  pattern,  and  another  with  the  fields  arranged 
in  a  V-shaped  pattern  (Fig.  2).  During  this,  as  well  as  any  other  Triangle/V-shaped  irradiation 
pair,  only  the  order  in  which  the  radiation  fields  are  delivered  changes.  Following  irradiation, 
cells  were  incubated  for  10-14  days,  stained,  and  counted  as  described  above  in  Subtask  LA. 
Surviving  fraction  was  also  calculated  as  described  above. 

Analysis  of  the  parameter  space  of  the  modeling  study  was  achieved  through  three 
separate  but  complementary  experiments.  The  high  a/|3  SQ-20  cells  were  included  in  the 
experiments  just  described  to  determine  the  effects  of  changing  tissue  type.  Then,  separate  plates 
of  WiDr  cells,  the  lowest  cx/|3  and  highest  P  line,  were  subjected  to  two  additional  radiation 
patterns.  In  one  of  these,  two  WiDr  plates  were  irradiated  with  the  same  patterns  and  Tf  as  in 
Fig.  1  with  each  of  the  six  fields  reduced  by  80%,  resulting  in  a  total  dose  of  180  cGy;  by 
comparing  this  to  the  900  cGy  results,  the  effect  of  dose/fx  was  analyzed.  Then  two  plates  were 
irradiated  with  the  same  patterns  and  doses/field  as  in  Fig.  1,  but  with  Tf  =5  min,  measuring  the 
effect  of  changing  Tf. 

The  analysis  described  in  Subtask  LB  Objective  2  was  done  statistically  on  each 
individual  experiment.  For  each  set  of  irradiation  conditions  (cell  line,  Tf,  dose/fx),  the  surviving 
fractions  of  cells  from  each  well  in  the  plate  irradiated  with  the  triangle  irradiation  pattern  was 
compared  to  the  measured  surviving  fractions  from  the  corresponding  V-shaped  pattern 
irradiated  plate  through  p-values  calculated  with  the  one-tailed  Mann-Whitney  U-test  [35]. 
Statistically  significant  differences  were  assumed  at  the  p<0.05  level. 

An  average  normalized  surviving  fraction  for  irradiation  with  the  Triangle  (Sto)  and  V- 
shaped  (Sv)  patterns  was  calculated  for  each  individual  experiment  by  averaging  the  measured 
normalized  surviving  fraction  from  all  six  wells  in  the  corresponding  plates  for  that  experiment. 
The  standard  error  in  each  Sxri  and  Sv  was  calculated,  using  standard  error  propagation 
techniques,  from  the  standard  deviation  of  the  PE  value  for  that  experiment  and  the  standard 
deviation  of  the  corresponding  non-normalized  surviving  fractions  used  to  calculate  each  factor. 
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Figure  2:  a)  A  schematic  of  the  six  field  900  cGy  dose  pattern  used  to  irradiate 
cells.  Here  the  fields  have  been  arranged  in  the  “Triangle”  pattern  (dashed  line) 
which  yields  maximum  cell  kill.  In  the  figure,  the  horizontal  axis  represents  time 
b)  The  same  as  (a),  except  a  representation  of  the  “V-shaped”  dose  pattern 
(dashed  line)  to  yield  minimum  cell  kill.  Note  the  exact  same  fields  are  used  in 
both  (a)  and  (b)  and  the  duration  of  delivery  is  exactly  the  same  -  only  the  order 
in  which  they  were  delivered  has  been  changed. 
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All  three  cell  lines  analyzed  were  irradiated  with  the  dose  patterns  in  Fig.  2  to  a  total  dose 
of  900  cGy  over  20  min.  Figs.  3a-c  show  the  measured  Stu  and  Sv  values  for  PC-3,  WiDr,  and 
SQ-20B,  respectively,  for  all  three  experimental  iterations  on  each  cell  line.  Above  each  Sxri/Sv 
pair  is  the  corresponding  p-value  for  that  experiment.  Note  that  for  the  high-|3  cell  lines,  WiDr 
(Fig.  3a)  and  PC-3  (Fig.  3b),  there  is  a  consistent,  visible  increase  in  surviving  fraction  when 
irradiating  with  the  V-shaped  pattern  versus  the  Triangle  pattern.  Furthermore,  p<0.05  for  each 
individual  WiDr  and  PC-3  experiment,  indicating  that  the  increase  in  survival  from  Triangle  to 
V-shaped  irradiation  is  statistically  significant  in  all  cases.  These  results  show  that  tangible 
differences  in  cell  survival  can  be  achieved  simply  by  rearranging  the  temporal  pattern  of  applied 
dose. 
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Figure  3:  Comparison  of  surviving  fraction  of  cells  following  irradiation  with  the  Triangle  dose  pattern  versus  the 
V-shaped  dose  patterns  as  discussed  in  the  text  using  a  total  dose  900  cGy  and  fraction  duration  of  20  min.  a)  WiDr 
cells,  b)  PC-3  cells,  and  c)  SQ-20B  cells.  For  each  combination  of  cell  line  and  irradiation  condition,  three 
independent  experimental  iterations  were  performed  (marked  Experiments  1-3  on  graphs).  One  tailed  p-values 
calculated  using  the  Mann-Whitney  U-test  for  all  experimental  iterations  in  all  cell  line-irradiation  condition 
combinations  are  shown  above 
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The  high-cx/p  SQ-20B  cells  show  a  radically  different  picture.  In  Fig.  3c  the  visible 
differences  seen  for  the  WiDr  and  PC-3  cells  are  minimal  or  nonexistent.  The  comparative 
statistics  for  each  SQ-20B  experiment  either  did  not  show  statistically  significant  differences 
between  the  populations  (p>0.05  for  iterations  1  and  3)  or  borderline  statistically  significant 
differences  (p~0.05  for  iteration  2).  This  loss  of  persistent,  statistically  significant  differences  for 
the  same  irradiation  patterns  and  parameters  applied  to  the  SQ-20B  cells  agrees  with  the 
assertion  of  the  modeling  study  that  differences  in  cell  survival  between  irradiation  with  the 
Triangle  and  V-shaped  patterns  decrease  with  increasing  cx/p  [34]. 

Figs.  4a, b  show  the  same  comparisons  for  the  WiDr  cells  when  irradiated  with  180  cGy 
over  20  min,  or  900  cGy  over  5  min,  respectively.  The  consistent  increase  in  survival  when 
irradiating  with  the  V-shaped  dose  pattern  compared  to  the  Triangle  dose  pattern  as  seen  for  the 
900  cGy/20  min  irradiation  in  Fig.  3  a  is  absent.  Furthermore  for  all  experimental  iterations  of 
both  the  180  cGy/20  min  and  900  cGy/5  min  irradiations  the  populations  receiving  Triangle  and 
V-shaped  dose  patterns  do  not  show  statistically  significant  differences  (p>0.05).  These  results 
both  agree  with  the  modeling  study,  which  predicted  that  the  increase  in  cell  survival  between 
irradiation  with  the  Triangle  versus  V-shaped  dose  patterns  is  larger  with  greater  dose/fx  and/or 
Tf[34]. 
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Figure  4:  Similar  comparison  to  Fig.  3a  of  surviving  fraction  of  WiDr  cells  following  irradiation  with  the  Triangle 
dose  pattern  versus  the  V-shaped  dose  patterns  as  discussed  in  the  text  for  a)  a  total  dose  (D)  of  180  cGy  and  fraction 
duration  (Tf)  of  20  min  b)  D  =900  cGy  and  Tf  =5  min.  Unlike  Fig.  3a,  p-values  above  each  corresponding 
Triangle/V-shaped  survival  pair  do  not  show  statistically  significant  differences. 
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We  then  eompleted  Subtask  LB  Objective  2  by  calculating  the  average  difference  in 
surviving  fraction  when  irradiating  with  the  Triangle  versus  V-shaped  dose  pattern  for  each  set  of 
irradiation  conditions.  These  were  then  compared  to  corresponding  values  predicted  by  the 
modeling  study.  The  average  measured  percent  change  in  survival  between  irradiation  with  the 
Triangle  and  V-shaped  pattern  (A%STri-v)  was  calculated  for  each  set  of  irradiation  conditions  by: 

•  100  (4) 

n=\  ^V,n 

where  n  =  1,2,3  represent  the  three  experimental  iterations  performed  for  each  set  of  conditions. 
Standard  errors  in  A%STri-v  were  again  calculated  using  standard  propagation  methods. 

This  data  was  compared  to  that  predicted  by  the  theoretical  models  by  substituting  Eq.  1 
into  Eq.  3  (ignoring  the  summation): 

=  [l  -  exp(;^(G^„.  (t)  -  Gy  (t))D^  )J- 100  (5) 

where  GTri(t)  and  Gv(t)  were  calculated  using  the  corresponding  temporal  patterns,  Eq.  2,  and  the 
techniques  (including  population-average  tissue-dependent  parameters)  of  the  modeling  study 
[34].  We  used  the  population  average  parameters  to  calculate  the  G(t)  values  for  reasons 
previously  discussed  in  Subtask  LA.  The  (3  values  were  those  we  previously  measured. 

These  relationships  between  the  measured  and  calculated  A%STri-v  values  can  be  seen  in 
Eig.  5.  A%STri-v  was  21.2-i-/-2.5%  for  the  WiDr  cells  and  18.6-i-/-4.7%  for  the  PC-3  cells  at  900 
cGy/fx  and  Tf=20  min.  Eor  the  higher  cx/p  SQ-20B  cells  this  drops  to  A%STri-v=3.0-i-/-1.3%. 
Irradiating  WiDr  cells  with  180  cGy/fx  over  20  min  resulted  in  A%STri-v=2.7-i-/-4.1%,  while 
irradiating  them  with  900  cGy/fx  over  5  min  yielded  A%STri-v=-0.8-i-/-4.2%.  Eig.  5  shows  that 
these  results  compare  very  favorably  to  those  predicted  by  the  model,  although  the  model- 
predicted  values  should  only  be  seen  as  rough  estimates  due  to  the  errors  in  the  measured  P 
values  and  the  use  of  population  average  parameters  in  the  calculation  of  GTri(t)-Gv(t)  (Eq.  4). 
Ultimately,  measurement  and  theory  display  the  same  patterns  throughout:  relatively  large 
increase  when  comparing  the  surviving  fraction  of  cells  irradiated  with  the  Triangle  dose  pattern 
versus  the  V-shaped  for  low-a/p  cell  lines,  large  doses/fx  and  long  Tf.  However,  if  oc/p  is 
decreased,  the  dose/fx  is  decreased,  or  Tf  is  shortened,  these  differences  diminish. 
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Figure  5:  Average  percent  change  in  surviving  fraction  (%AS)  between  irradiation  with  the  Triangle  and  V-Shaped 
dose  pattern  for  each  cell  line  and  irradiation  condition  combination  shown  in  Figs.  4  and  5  (blue  bar).  Error  bars  on 
these  curves  indicate  the  standard  error  in  the  %AS  amongst  the  three  experimental  iterations  performed  for  each  set 
of  conditions.  Yellow  bars  indicate  the  %AS  values  predicted  using  the  measured  oc/p  values  for  each  cell  line  (see 
text)  along  with  other  average  population  parameters  [34]. 


Task  11:  Characterize  a  phantom  for  use  in  three-dimensional  in  vitro  experiments  in  a  realistic 

treatment  environment 


II.  A:  Characterization  of  dosimetry  in  IMRT  radiobiological  experiment  phantom  using  TLDs 
and  film.  (7-10  mos.) 

Objectives: 

1.  Create  IMRT  treatment  plan(s)  for  phantom.  (7-8  mos.) 

2.  Determine  consistently  responding  TLD  set  for  experiment.  (7-8  mos.) 

3.  Characterize  phantom  dosimetry  with  TLDs  and  film.  (8-10  mos.) 

4.  Analysis  of  measured  dosimetry  with  TLDs  and  film  compared  to  predicted  dosimetry 

from  treatment  planning  software.  Determine  a  scheme  for  reconciliation  of  the 
two  as  needed.  (9-10  mos.) 

All  of  the  objectives  of  subtask  II.A  were  completed.  The  data  and  analysis  were  published  as 
Altman  et  al.  ’’Characterization  of  a  novel  phantom  for  three-dimensional  in  vitro  cell 
experiments.”  Phys  Med  Biol.  2009;  54(5):N75-N82.  This  paper  has  been  attached  in  the 
appendix  of  this  report;  it  will  provide  an  in  depth  summary  and  analysis  of  the  completion  of  the 
objectives  in  subtask  II.A. 

Il.B:  Assess  the  effectiveness  of  the  phantom  as  a  tool  for  performing  in  vitro  experiments. 

(10-14  mos.) 

Objectives: 
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1  Experimental  design  and  setup  (10-11  mos.) 

2.  Perform  in  vitro  experiments,  repeating  experiments  to  ensure  result  eonsisteney  as 

needed.  (11-13  mos.) 

3.  Analyze  results  to  determine  if  equivaleney  with  standard  in  vitro  experimental 

teehniques  is  feasible  (13-14  mos.) 

The  effeetiveness  of  the  phantom  as  a  tool  for  performing  in  vitro  experiments  was 
assessed  through  a  biological  validation  using  two  human  cancer  A549  [36,  37],  a  lung 
adenocarcinoma  (obtained  from  the  American  Type  Culture  Collection),  and  SCC116  [38],  an 
alveolar  ridge  squamous  cell  carcinoma  (acquired  from  the  lab  of  Dr.  Susanne  Gollin  at  the 
University  of  Pittsburgh).  Cells  were  analyzed  using  the  diphenylamine  (DPA)  assay,  which 
measures  DNA  accumulation  or  loss,  and  can  be  used  to  calculate  the  percent  of  viable  cells 
post-irradiation  (%Via)  [39].  In  contrast  to  the  clonogenic  survival  assay,  the  colorimetric  nature 
of  the  DPA  assay  means  that  more  data  points  can  be  efficiently  collected.  Furthermore,  the 
automated  readout  systems  used  to  analyze  DPA  assays  limit  the  human  error  often  associated 
with  counting  colonies  in  the  clonogenic  assay.  All  chemicals  used  in  the  preparation  of 
solutions  needed  for  the  assay  were  purchased  from  the  Sigma-Aldrich  Chemical  Company  and 
used  as  received.  Cells  were  seeded  (100  pL)  into  96- well  plates  and  grown  until  they  were 
approximately  50-60%  confluent.  Nine  identical  plates  were  prepared,  with  one  (randomly 
selected)  of  these  serving  as  the  unirradiated  control.  Three  hours  prior  to  irradiation,  an 
additional  150  pL  of  media  was  added  to  each  well. 

The  plates  were  irradiated  using  two  different  experimental  setups:  a  “standard”  setup 
similar  to  that  used  in  a  typical  radiobiological  experiment  described  earlier,  and  one  employing 
the  IMRT  phantom.  In  both  cases,  doses  of  I,  2,  5,  and  10  Gy  were  delivered  to  the  cell  lines 
using  the  linear  accelerator.  In  the  standard  setup,  each  dose  was  imparted  using  a  single  gantry 
angle  (0°)  to  one  96-well  plate  placed  within  a  solid  water  block.  6  MV  photons  were  delivered 
at  600  MU/min  with  the  well  plate  at  the  center  of  a  20  x  20  cm^  field  at  90  cm  SSD  and  10  cm 
depth.  For  the  IMRT  setup,  a  treatment  plan  was  created  which  contained  two  -235  cm  PTVs, 
each  containing  four  rows  of  eight  wells  of  all  plates.  Six  equally-spaced  beam  angles  were  used, 
ranging  from  28°-336°.  Each  PTV  resided  at  opposite  ends  of  the  well  plates.  The  ratio  of  doses 
between  the  two  PTVs  was  fixed  at  2:1.  Two  plans  were  created  to  deliver  nominal  doses  of  1 
and  2  Gy,  and  5  and  10  Gy  to  the  proximal  and  distal  PTVs,  respectively.  The  phantom  was 
loaded  with  two  well  plates  per  plan. 

Immediately  following  irradiation,  the  plates  were  placed  back  in  the  incubator  for  an 
additional  96  hours.  After  this  time,  the  media  was  completely  removed  from  the  cells,  and  60  |a,F 
of  a  1:5  mixture  of  chilled  acetaldehyde  (0.16%  in  water)  and  perchloric  acid  (20%  vol/vol)  was 
added.  Each  well  was  then  treated  with  100  jiL  of  a  4%  DPA  solution  in  glacial  acetic  acid,  and 
the  plates  were  incubated  for  an  additional  24  hours  at  37  °C.  The  absorbance  was  then  read  at 
595  nm  for  each  well.  Individual  trials  were  normalized  to  mean  readings  from  the  unirradiated 
control  plates  to  calculate  the  %Via  values,  which  were  then  averaged. 

For  each  cell  line,  single  iterations  of  the  biological  validation  experiment  were  analyzed 
separately.  The  error  in  every  measured  %Via  was  determined  from  the  average  a  of  all  the 
points  under  the  same  irradiation  conditions  and  that  of  the  unirradiated  control  readings  using 
standard  error  propagation  techniques.  The  average  %Via  standard  deviation  for  each  irradiation 
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condition  of  the  experiment  was  then  ealeulated.  For  a  given  dose  point,  the  %Vias  generated 
using  the  two  different  experimental  setups  were  compared  using  student  t- tests. 

For  both  the  IMRT  phantom  and  the  standard  setup,  average  %Vias  and  assoeiated 
standard  deviations  from  three  independent  experiments  was  ealeulated  for  both  SCC116  and 
A549  eells  at  several  dose  points  from  0-10  Gy  (Table  1).  The  differences  in  %Via  averaged 
1.28%  and  3.26%  for  the  SCC116  and  A549  eells,  respeetively.  The  results  for  the  two 
experimental  setups  revealed  p>0.05  when  eompared  with  the  student  t-test  for  both  eell  lines  at 
all  dose  points.  As  the  differenees  in  %Via  do  not  show  statistieal  signifieanee  for  either  eell 
line,  this  implies  that  experimental  results  acquired  in  the  phantom  eould  be  eompared  to  or 
evaluated  against  results  from  previous  or  future  studies  employing  more  standard  teehniques. 


Table  1:  Percentage  of  viable  cells  post-irradiation  measured  using  the  DPA  assay 
for  two  different  cell  lines. 


Percentage  of  Viable  Cells  Post-Irradiation  (%) 

Cell 

Dose 

Standard  Setup 

IMRT  Phantom 

Absolute  ^ 
Difference 

Line 

(cGy) 

Average 

Standard 

Error 

Average 

Standard 

Error 

100 

89.23 

8.04 

87.40 

8.91 

1.84 

SCC116 

200 

72.25 

9.45 

73.15 

7.38 

0.90 

500 

31.15 

5.51 

30.38 

3.71 

0.76 

1000 

12.30 

3.28 

13.94 

4.20 

1.63 

100 

55.35 

8.66 

60.59 

12.38 

5.24 

A549 

200 

34.49 

7.74 

38.22 

8.15 

3.73 

500 

14.24 

4.24 

17.59 

5.14 

3.35 

1000 

13.64 

8.76 

12.93 

5.04 

0.71 

Abbreviations:  DPA  =  diphenylamine  IMRT  =  Intensity  modulated  radiation  therapy 
*  Absolute  difference  is  between  the  average  percentage  of  viable  cells  post-irradiation  from  the  standard  setup  and 
IMRT  phantom.  At  each  dose  point,  all  differences  show 
p  >  0.05. 


KEY  RESEARCH  ACCOMPLISHMENTS 


•  Determined  cell  lines  of  interest  for  in  vitro  experiments. 

•  Showed  the  effect  of  the  rearrangement  of  the  temporal  pattern  of  applied  dose  for 
a  single  fraction  of  radiation  in  vitro  is  statistically  significant  for  low  oc/p  cell 
lines  at  high  dose/fx  and  long  fraction  duration,  as  predicted  by  the  previously 
published  modeling  study. 
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•  Showed  that  temporal  pattern  based  statistically  significant  differences  disappear 
in  vitro  with  increasing  cx/|3,  decreasing  dose/fx,  or  decreasing  fraction  duration, 
results  which  also  agree  with  conclusions  of  the  modeling  study. 

•  Verified  dosimetrically  that  patient-based  treatment  plans  can  be  accurately 
projected  into  a  new  phantom  for  use  in  three  dimensional  in  vitro  experiments. 

•  Performed  a  biological  validation  study  showing  that  the  use  of  the  phantom  is 
equivalent  to  established  radiobiological  experimental  techniques. 


REPORTABLE  OUTCOMES 


Original  Peer-reviewed  Publications: 

Altman  MB,  Vesper  BJ,  Smith  BD,  Stinauer  MA,  Pelizzari  CA,  Aydogan  B,  Reft  CS, 
Radosevich  JA,  Chmura  SJ,  Roeske  JC.  Characterization  of  a  novel  phantom  for  three- 
dimensional  in  vitro  cell  experiments.  Phys  Med  Biol.  2009;  54(5):N75-N82. 

Altman  MB,  Stinauer,  MA,  Javier,  D,  Aydogan  B,  Pelizzari,  C,  Chmura  SJ,  Roese,  JC. 
Validation  of  Temporal  Optimization  for  a  Single  Traction  of  Radiation  in  Vitro.  Int  J  Rad  One 
Biol  Phys  (Submitted). 

Abstracts  Presented  at  National  or  International  Scientific  Conferences: 

Altman  MB,  Vesper  BJ,  Smith  BD,  Stinauer,  MA,  Pelizzari  CA,  Aydogan  B,  Radosevich  JA, 
Chmura  SJ,  Roeske  JC.  Dosimetric  Characterization  and  Biological  Validation  of  a  Novel 
Phantom  for  Three-Dimensional  IMRT-Based  in  Vitro  Experiments.  Med  Phys.  35(6):  2816. 
2008 

Hoggarth  MA,  Stinauer  MA,  Altman  MB,  Roeske  JC.  Automated  Colony  Counting  Using  Color 
and  Image  Processing  Techniques.  Med  Phys.  35(6):  2816.  2008 

CONCLUSION 


Through  a  series  of  in  vitro  experiments,  we  have  verified  several  predictions  of  a 
previously  performed  modeling  study.  We  have  demonstrated  that,  for  a  single  fraction  of  dose, 
cell  survival  can  be  affected  simply  by  manipulation  of  the  temporal  pattern  of  dose  applied 
during  that  fraction.  It  was  then  shown  that  these  differences  are  greatest  in  situations  with  low- 
cx/p  and/or  high-P  tissues,  high  doses/fx  and  long  treatment  times.  Although  the  effects  of 
temporal  dose  manipulation  need  to  be  further  analyzed  in  vivo,  the  results  of  both  the  present 
work  and  the  modeling  study  provide  a  template  of  how  temporal  optimization  techniques  could 
potentially  be  used  to  positively  aid  treatment  outcome,  either  through  increasing  cell  kill  to 
tumor  tissue  or  enhancing  normal  tissue  sparing. 
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To  test  temporal  optimization  effeets  in  a  realistie  treatment  environment  we  have  created 
a  simple,  flexible,  and  efficient  phantom  for  use  in  three-dimensional  IMRT-based  in  vitro 
experiments.  The  simultaneous  use  of  multiple  tissue-culture  plates  provides  a  three- 
dimensional  grid  of  biological  sample  points.  By  projecting  and  delivering  three-dimensional 
dose  distributions  into  the  cylindrical  phantom,  radiobiological  experiments  can  be  performed  in 
a  realistic  IMRT  setting.  Multiple  points  within  the  dose  delivery  volume  can  be  analyzed  in  a 
single  irradiation  session.  The  TLD  and  film  analysis  performed  in  this  study  show  that  IMRT- 
based  dose  distributions  can  be  accurately  and  reliably  projected  into  the  phantom  under  a  wide 
range  of  experimental  designs.  Equivalent  biological  outcome  over  a  range  of  doses  for  cells 
irradiated  within  the  phantom  versus  those  irradiated  in  a  more  standard  radiobiological  setup 
show  that  the  phantom  is  a  logical  progression  in  the  experimental  investigation  of 
radiobiological  effects  in  emerging  radiation  dose  delivery  technologies. 

With  the  movement  towards  patient-specified  treatment  for  prostate  cancer,  including 
biological  parameters  and  influences,  the  results  of  the  in  vitro  study  showed  that  temporal 
optimization  techniques  could  prove  to  be  an  important  element  in  enhancing  the  efficacy  of 
radiation  therapy  for  prostate  cancer.  The  dosimetric  accuracy  and  experimental  effectiveness  of 
the  phantom  make  it  a  good  potential  tool  for  analyzing  temporal  optimization  effects  in  vitro 
using  realistic  clinical  prostate  cancer-based  treatments. 
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Abstract 

A  novel  intensity-modulated  radiation  therapy  (IMRT)  phantom  for  use  in 
three-dimensional  in  vitro  cell  experiments,  based  on  a  commercially  available 
system  (CIRS  Inc.,  Norfolk,  VA),  was  designed  and  fabricated.  The  water- 
equivalent  plastic  phantom  can,  with  a  set  of  water-equivalent  plastic  inserts, 
enclose  1-3  multi-well  tissue  culture  plates.  Dosimetry  within  the  phantom  was 
assessed  using  thermoluminescence  dosimeters  (TLDs)  and  film.  The  phantom 
was  loaded  with  three  tissue  culture  plates,  and  an  array  of  TLDs  or  a  set  of  three 
films  was  placed  underneath  each  plate  within  the  phantom,  and  then  irradiated 
using  an  IMRT  plan  created  for  it.  Measured  doses  from  each  dosimeter  were 
compared  to  those  acquired  from  the  treatment  planning  system.  The  percent 
differences  between  TLD  measurements  and  the  corresponding  points  in  the 
treatment  plan  ranged  from  1.3%  to  2.9%,  differences  which  did  not  show 
statistical  significance.  Average  point-by-point  percent  dose  differences  for 
each  film  plane  ranged  from  1.6%  to  3.1%.  The  percentage  dose  difference 
for  which  95%  of  the  points  in  the  film  matched  those  corresponding  to  the 
calculated  dose  plane  to  within  3.0%  ranged  from  2.8%  to  4.2%.  The  good 
agreement  between  predicted  and  measured  dose  shows  that  the  phantom  is  a 
useful  and  efficient  tool  for  three-dimensional  in  vitro  cell  experiments. 


1.  Introduction 

Intensity  modulated  radiation  therapy  (IMRT)  is  a  widely  used  form  of  external  beam  radiation 
therapy  (EBRT)  in  which  the  beams  are  modulated  during  the  course  of  therapy.  The 
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modulation  provides  the  capability  to  produce  radiation  fields  of  non-uniform  intensity.  IMRT 
allows  the  delivered  dose  distributions  to  conform  more  closely  to  the  volume  of  the  target  than 
could  generally  be  achieved  through  conventional  EBRT.  The  volumetric  conformity  attainable 
with  IMRT  results  in  a  highly  variable  dose  distribution,  both  spatially  and  temporally.  An 
increasing  number  of  studies  have  investigated  the  effects  of  applying  dose  using  a  variable 
spatial  and/or  temporal  dose  pattern  on  the  underlying  biology  in  vitro  (Mu  et  al  2003, 
Suchowerska  et  al  2005,  Bromley  et  al  2006,  Mackonis  et  al  2007,  Sterzing  et  al  2005, 
Moiseenko  et  al  2007). 

However,  testing  the  effects  of  an  IMRT-style  dose  distribution  in  a  biological  setting 
is  not  straightforward.  Many  experiments  use  a  setup  where  a  single  cell-containing  vessel 
is  encased  in  a  block  of  water-equivalent  material  and  irradiated  using  a  single  fixed  beam 
(Mu  et  al  2003,  Suchowerska  et  al  2005,  Bromley  et  al  2006,  Mackonis  et  al  2007).  The  dose 
distributions  in  these  setups  are  essentially  two  dimensional,  in  that  the  population  of  cells  is 
irradiated  with  a  single  dose  pattern.  However,  in  an  IMRT-based  experiment,  there  is  a  need 
to  analyze  cell  response  at  discrete  positions  within  the  entire  three-dimensional  treatment 
volume.  In  one  study,  Sterzing  et  al  (2005)  created  a  specialized  cylindrical  phantom  to  try 
and  resolve  these  issues.  It  used  an  array  of  nine  co-planar  sample  points  acquired  using 
cryotubes;  the  position  of  the  sample  plane  was  translatable  both  in  the  axial  and  lateral 
dimensions  allowing  for  sampling  of  a  full  three-dimensional  space.  However,  the  radiation 
plan  would  need  to  be  applied  at  every  position  of  the  array  necessary  to  sample  the  dose 
space,  limiting  its  experimental  efficiency.  Furthermore,  irradiation  in  cryotubes  requires 
plating  cells  into  tissue  culture  plates  or  flasks  before  data  acquisition.  This  introduces  more 
potential  experimental  error  into  the  results,  due  in  part  to  the  subjugation  of  the  cells  to  an 
additional  environmental  change. 

An  optimal  tool  would  (1)  provide  the  capability  to  perform  in  vitro  cell  experiments 
three-dimensionally  in  an  IMRT  setting,  (2)  allow  for  efficient  data  collection  in  all  three 
dimensions  of  the  treatment  volume,  and  (3)  minimize  experimental  error  by  limiting  the 
number  of  necessary  cell  transfers.  We  have  created  a  technique  using  a  simple  custom- 
designed  phantom  which  meets  all  three  of  the  above  challenges. 

2.  Methods  and  materials 

2.7.  IMRT  phantom 

A  cylindrical  phantom  was  designed  which  is  30  cm  long  with  a  16  cm  diameter,  based  on  a 
commercially  available  neck  phantom  (Model  002HN,  CIRS  Inc.,  Norfolk,  VA).  The  phantom 
is  composed  of  a  proprietary  water-equivalent  plastic  (CIRS  Inc.,  Norfolk,  VA)  and  contains  a 
custom-cut  86  x  66  mm^  rectangular  bore  extending  half  the  phantom’s  length.  A  set  of  128  x 
86  mm^  rectangular  inserts,  ranging  from  1  to  20  mm  thick,  and  an  86  x  66  x  35  mm^  cap, 
all  made  from  the  same  water-equivalent  plastic  as  the  rest  of  the  phantom,  allow  up  to  three 
stacked  standard  multi-well  tissue  culture  plates  to  be  sealed  inside  the  phantom  with  minimal 
air  gaps  (figure  1(a)).  The  number  of  wells  within  a  plate  can  vary,  as  can  plate  thickness 
between  manufacturers,  although  a  rectangular  128  x  86  mm^  base  is  standard.  As  a  result,  the 
inserts  provide  the  phantom  with  a  flexibility  to  manage  a  wide  range  of  experimental  setups 
(figure  1(b)).  To  ensure  spatial  precision,  a  series  of  six  external  guide  markers  at  two  axial 
positions  on  the  phantom  are  aligned  to  the  in-room  laser  system.  Movement  during  scanning 
and  treatment  is  prevented  by  placing  the  phantom  on  a  plastic  stand  molded  to  the  curvature 
of  the  phantom.  Both  the  stand  and  the  guide  markers  are  components  of  the  original  CIRS 
Model  002HN  phantom. 
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Figure  1.  (a)  Phantom  schematic.  Multi-well  cell  culture  plates  can  he  placed  into  the  cavity 
during  experiments.  The  zone  labeled  ‘back  half  contains  markers  and  chambers  for  use  with 
an  array  of  dosimeters  and  is  the  same  as  that  in  the  Model  002HN  CIRS  phantom,  (b)  Phantom 
loaded  with  tissue  culture  plates  and  inserts,  (c)  Single  six-well  tissue  culture  plate  and  lid  (front). 
Stack  of  three  six-well  plates  (back). 

(This  figure  is  in  colour  only  in  the  electronic  version) 


2.2.  Dosimetric  characterization 

The  phantom  is  designed  to  work  with  standard  multi-well  tissue  culture  plates,  which  are 
generally  single  pieces  of  molded  polystyrene  plastic  with  separate  fitted  lids  (figure  1(c)). 
The  wells  are  a  rectangular  array  of,  essentially,  ‘mini  Petri  dishes’  recessed  from  the  top  of 
the  plate.  Under  experimental  conditions,  the  stacked  well  plates  present  a  matrix  of  plastic, 
liquid  (filling  the  wells),  and  air  pockets.  Hence,  the  factor  which  could  potentially  have 
the  greatest  adverse  affect  on  accurate  dose  delivery  within  the  phantom  under  experimental 
conditions  is  these  air  gaps  within  the  plates  between  the  cells  and  the  phantom  walls.  To 
assess  the  maximum  impact  of  these  air  gaps,  phantom  dosimetry  was  characterized  with  both 
thermoluminescence  dosimeters  (TLDs)  and  film  using  a  ‘worst  case’  experimental  setup  in 
terms  of  air  in  the  plate  containing  chamber.  A  clonogenic  survival  assay  setup,  resulting  in 
air  occupying  ~75%  of  the  plate  chamber  volume,  was  used;  a  stack  of  three  2.2  cm  thick 
six-well  tissue  culture  plates  (Coming  Inc.,  Corning,  NY)  filled  with  2  mL  H2O  per  well.  A 
CT  (PQ5000,  Philips  Medical  Systems,  Andover,  MA)  scan  of  the  plate-filled  phantom  was 
obtained.  An  IMRT  plan  for  this  CT  scan  was  generated  using  the  Eclipse  treatment  planning 


N78 


M  B  Altman  et  al 


6  Well  Plate  To  View 


TLD  Pair 

I  I  =  Plates  1 , 2,  and  3 
:  -I  =  Plate  2  Only 


Phantom  Axial  View 


Figure  2.  Schematic  six-well  plate  placement  within  the  phantom  (left)  and  placement  of  ten  TLD 
pairs  on  the  bottom  of  plates  I,  2  and  3  within  the  phantom  (right).  Dark  grey  rectangles  denote 
TLD  locations  common  to  all  three  plates.  Dotted  rectangles  denote  TLD  locations  on  plate  2  only. 


system  (Varian  Corporation,  Palo  Alto,  CA)  containing  a  single  planning  target  volume  (PTV) 
encompassing  all  three  plates  with  mean  dose  200.9  cGy  and  —  4.53  cGy.  Six  equally 
spaced  beam  angles  were  used,  ranging  from  28°  to  336°.  Dose  calculations  were  performed 
using  the  analytical  anisotropic  algorithm  (AAA)  (Ulmer  and  Kaissl  2003,  Ulmer  et  al  2005). 
All  irradiations  in  this  study  were  performed  using  a  clinical  linear  accelerator  (2100CD, 
Varian  Medical  Systems,  Palo  Alto,  CA). 

A  set  of  TLDs  with  dimensions  0.038  x  0.32  x  0.32  cm^  and  sensitivities  with  standard 
deviations  <1%  was  determined  using  the  methods  described  by  Al  Hallaq  et  al  (2006)  with 
a  Harshaw  3500  TLD  Reader  (Thermo  Fisher  Scientific  Inc.,  Waltham,  MA).  Individual  TLD 
sensitivity  correction  factors  were  determined  from  this  process  and  used  to  correct  all  the 
subsequent  TLD  measurements.  Pairs  of  TLDs  were  placed  at  ten  locations  among  the  three 
water-filled  plates,  which  were  then  loaded  in  the  phantom  (figure  2).  The  TLDs  were  placed 
on  the  underside  of  the  water-filled  wells  in  the  plate  as  cells  adhere  to  the  floor  of  the  wells 
in  the  clonogenic  assay.  Such  a  placement  puts  the  TLDs  (as  with  the  film,  discussed  below) 
as  close  to  the  position  of  the  cells  under  experimental  conditions  as  possible  while  keeping 
them  dry. 

The  phantom  was  irradiated  with  the  clinical  linear  accelerator  using  the  IMRT  plan.  This 
was  repeated  twice,  each  time  placing  the  same  TLD  at  the  same  position  in  the  phantom. 
Each  IMRT  delivery  included  calibrated  TLDs  to  convert  the  TLD  readings  to  absorbed 
dose.  The  predicted  doses  to  the  TLDs  were  calculated  in  Eclipse.  Systematic  errors  in 
the  TLD  measurements  were  determined  to  be  1.5%  from  the  initial  calibration.  Individual 
TLD  random  errors  were  assumed  to  be  the  standard  deviations  of  the  three  independent 
in-phantom  measured  doses.  The  systematic  error  in  the  treatment  planning  calculated  doses 
was  estimated  at  3.0%  (Khan  2003,  Aydogan  et  al  2007,  Sanchez-Doblado  et  al  2007).  The 
random  errors  in  calculated  dose  were  determined  from  standard  deviations  of  dose  in  the 
volumes  used  to  determine  the  expected  TLD  doses.  The  total  errors  were  determined  by 
adding  the  corresponding  random  and  systematic  errors  in  quadrature.  The  Student’s  f-test 
was  then  used  to  analyze  the  difference  between  measured  and  calculated  TLD  doses. 
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Table  1.  Comparison  between  measured  and  calculated  TLDs  dose. 


TLD  number 

Measured  dose  (cGy) 

Calculated  dose  (cGy) 

Percent  difference‘s 

Average 

cr 

Average 

cr 

1 

206.2 

1.1 

208.2 

0.3 

1.0% 

2 

205.4 

1.0 

208.1 

0.3 

1.3% 

3 

200.3 

1.8 

203.8 

0.3 

1.7% 

4 

199.2 

2.0 

203.7 

0.2 

2.2% 

5 

197.8 

1.0 

203.6 

0.2 

2.9% 

6 

200.6 

1.3 

203.8 

0.2 

1.6% 

7 

204.1 

0.6 

205.8 

0.3 

0.8% 

8 

205.7 

1.0 

205.8 

0.3 

0.0% 

9 

200.2 

1.5 

202.7 

0.2 

1.2% 

10 

203.0 

0.8 

202.7 

0.1 

-0.2% 

11 

200.1 

0.2 

203.1 

0.4 

1.5% 

12 

199.5 

3.0 

202.7 

0.4 

1.6% 

13 

197.3 

1.4 

201.3 

0.2 

2.0% 

14 

198.3 

1.7 

201.4 

0.1 

1.5% 

15 

203.2 

3.0 

201.6 

0.4 

-0.8% 

16 

204.9 

1.0 

202.2 

0.4 

-1.3% 

17 

199.4 

4.0 

199.7 

0.4 

0.2% 

18 

200.2 

4.4 

200 

0.5 

-0.1% 

19 

199.3 

5.0 

197.4 

0.4 

-1.0% 

20 

197.8 

1.6 

197.9 

0.4 

0.0% 

Abbreviations;  TLD  =  thermoluminescent  dosimeter,  cr  =  standard  deviation. 

^  Percent  difference  is  between  average  measured  and  calculated  dose.  All  values  show  no 
statistically  significant  differences  {p  >  0.05). 


The  treatment  plan  was  then  delivered  with  Kodak  EDR2  film  (Eastman  Kodak  Company, 
Rochester,  NY)  cut  to  plate  dimensions  and  placed  below  each  well  plate,  yielding  three  films. 
Each  film  was  irradiated  individually  within  the  phantom,  as  a  1-3%  error  was  found  due  to 
the  large  volume  of  air  in  the  phantom  setup  when  the  films  were  irradiated  together  (data 
not  shown).  Eollowing  irradiation,  films  were  processed  using  a  Kodak  X-Omat  3000RA 
(Eastman  Kodak  Company,  Rochester,  NY)  and  digitized  with  an  Epson  Expression  lOOOOXL 
flat-bed  scanner  (Epson  America  Inc.,  Long  Beach,  CA).  The  phantom-irradiated  films  were 
then  registered  to  corresponding  dose  plane  images  extracted  from  Eclipse  using  pre-existing 
software  developed  at  the  University  of  Illinois  at  Chicago.  Due  to  the  lack  of  steep  dose 
gradients  within  the  treatment  plan,  the  absolute  percentage  dose  difference  (%Ddifd  between 
corresponding  pixels  in  the  film  and  predicted  dose  plane  images  was  used  to  compare  the 
films  in  this  study.  The  mean  %Ddiff  was  calculated  for  all  pixels  within  the  interior  90%  of 
each  film  image,  eliminating  errors  due  to  edge  effects.  Eor  each  film,  the  at  which  95% 

of  the  points  in  the  image  were  equal  to  or  below  that  cutoff  (%Ddiff  95)  was  also  calculated. 

3.  Results 

The  average  TLD  dose  from  three  independent  measurements  in  the  phantom  ranged  from 
197.3  to  206.2  cGy,  while  the  doses  at  those  points  predicted  by  the  treatment  planning  software 
ranged  from  197.4  to  208.1  cGy  (table  1).  The  signed  percent  differences  in  table  1  indicate 
whether  the  measured  TLD  dose  was  lower  (negative)  or  higher  (positive)  than  the  predicted 
dose  and  range  from  —1.33%  to  2.86%.  Eor  all  TLDs,  f-test  comparison  of  measured  and 
predicted  doses  revealed  p  >  0.05.  This  implies  that  the  aforementioned  percent  differences 
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Table  2.  Measured  doses  for  the  three  films  displayed  in  figure  2  and  calculated  doses  of  the 
corresponding  dose  planes  from  the  treatment  planning  system. 


Dose  (cGy) 

Measured  Calculated 

%£>d.ff(%) 

%Z)diff,95  (%) 

Mean 

a 

Mean 

a 

Mean 

a 

Film  1 

197.7 

2.3 

202.3 

3.9 

2.6 

0.8 

3.8 

Film  2 

195.7 

2.3 

201.5 

4.3 

3.1 

0.6 

4.2 

Film  3 

195.6 

2.5 

198.9 

3.7 

1.6 

0.7 

2.8 

Abbreviations.'  a  =  standard  deviation,  %/)diff  =  percentage  dose  difference,  %Ddiff,95  =  the 
percentage  dose  difference  at  which  95%  of  all  points  in  the  images  agree. 


do  not  show  statistical  significance.  Hence  the  doses  projected  into  the  phantom  match  those 
actually  delivered  at  the  TLD  locations  even  in  the  highly  inhomogeneous  environment  of  the 
well-plate-filled  phantom. 

Table  2  displays  the  average  dose  and  standard  deviation  of  dose  for  each  measured  film, 
those  same  factors  for  the  corresponding  treatment  plan  calculated  dose  planes,  as  well  as 
%Ddiff  comparisons  for  those  planes.  Note  the  good  general  agreement  for  all  three  films 
(%Z)diff  ranged  from  1.6%  to  3.1%,  while  ranged  from  2.8%  to  4.2%).  It  could  be 

argued  that  Film  3,  located  at  the  tissue  culture  plate/water-equivalent  plastic  (plate /plastic) 
interface  showed  the  better  agreement  than  the  other  two  films,  in  that  it  had  the  lowest 
measured  %Z)diff  and  %Z)diff,95  values.  However,  EDR2  film  can  be  an  inconsistent  dosimeter 
in  and  around  inhomogeneities  and/or  air  gaps  (Charland  et  al  2003,  Chetty  and  Charland 
2002,  Yoon  et  al  2007,  Gillis  et  al  2005,  Pai  et  al  2007).  The  use  of  prolonged  irradiation  with 
multiple  IMRT-type  fields  could  also  result  in  film  under  response  due  to  the  Schwarzschild 
effect  (Djouguela  et  al  2005).  Thus,  although  the  differences  increase  slightly  away  from  the 
plate /plastic  interface,  this  could  largely  be  explained  by  the  inaccuracies  of  film  as  a  dosimeter 
in  those  environments.  There  is  also  a  3.0%  error  in  the  treatment  planning  calculated  dose, 
as  noted  in  section  2.2.  Given  these  and  the  overall  small  magnitudes  of  %Z)diff  and  %Z)diff,95, 
the  measured  film  doses  agree  well  with  those  projected  into  the  phantom. 

4.  Discussion 

We  have  created  a  novel  phantom  to  easily  and  efficiently  perform  in  vitro  cell  experiments  in 
a  fully  three-dimensional  IMRT  setting.  Up  to  three  stacked  multi-well  tissue  culture  plates 
can  be  placed  inside  the  cylindrical,  water-equivalent  plastic  phantom  along  with  a  series 
of  inserts  to  minimize  unnecessary  air  gaps,  a  setup  providing  a  three-dimensional  array  of 
in  vitro  experiments.  A  strength  of  this  design  is  its  simplicity;  as  the  phantom  could  be 
easily  replicated  in  a  machine  shop,  many  research  centers  could  perform  experiments  using 
it.  Although  simple,  the  phantom  proves  advantageous  in  terms  of  experimental  technique.  As 
the  plate  space  is  large  compared  to  other  previous  experimental  designs,  the  phantom  provides 
a  far  more  complete  and/or  efficient  biological  sampling  of  a  three-dimensional  IMRT  dose 
space  than  any  of  the  previously  used  techniques.  Furthermore,  the  phantom  has  no  moving 
parts,  and  requires  no  specialized  equipment  or  water  filling  to  use,  enhancing  its  efficiency 
compared  to  other  methods. 

Measured  dosimetry  within  the  phantom  analyzed  under  experimental  conditions  using 
both  TLDs  and  film  agreed  well  with  the  dose  projected  into  the  phantom  by  the  treatment 
planning  system.  The  phantom  is  designed  for  potential  use  in  irradiation  schemes  with  large 
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dose  homogeneities.  We  did  not  irradiate  with  such  a  plan,  as  the  clonogenic  setup  we  have 
used  requires  large  sample  spaces  (i.e.  the  whole  well  would  be  irradiated  with  a  constant 
dose).  Furthermore,  the  configurations  with  sufficient  biological  sample  frequency  for  more 
complex  dose  distributions  would  use  either  plates  with  higher  number  wells /plate  or  ‘single 
weir  plates  (essentially  a  Petri  dish  with  multi-well  plate  base  dimensions).  These  plates  are 
generally  thinner  than  those  which  we  used  and,  when  filled,  would  result  in  far  less  air  by 
volume  in  the  plate  than  the  experimental  setup  we  employed.  This,  in  combination  with  the 
accuracy  in  the  worst  case  experimental  setup,  implies  the  phantom  could  be  reliably  used  in 
any  configuration  and  with  any  dose  distribution. 

As  well  plates  with  any  number  of  wells  will  fit  inside  the  phantom,  assays  which  require 
cells  to  adhere  to  the  well  floor  or  free  float  in  media  are  both  possible,  provided  a  well  of 
sufficient  size  is  chosen  to  allow  for  the  desired  dose  space  sampling.  All  of  these  factors 
mean  any  assay  performed  using  well  plates  is  feasible,  including  those  at  the  cutting  edge 
of  radiation  therapy  research.  Gene  microarry  analysis  of  genes  coding  for  radiosensitivity 
(Torres-Roca  et  al  2005,  Wong  et  al  2006),  radiation-based  studies  of  y— H2AX,  a  marker 
for  DNA  double  strand  breaks,  (Kao  et  al  2006)  and  several  enzyme-linked  immunosorbent 
assay  (ELISA)  studies  (Kaminski  et  al  2005,  Zhao  et  al  2008,  Fleckenstein  et  al  2007) 
have  all  addressed  a  number  of  radiation-therapy-related  issues.  The  phantom  provides  a 
simple,  efficient  platform  to  put  future  similar  experiments  into  a  realistic,  treatment-based 
environment;  patient-based  three-dimensional  dose  distributions  could  be  delivered,  allowing 
for  multiple  biological  samples  from  a  wide  array  of  points  within  the  treatment  volume. 

The  simplicity  and  flexibility  of  the  phantom  design  allow  it  to  be  used  with  a  wide  range 
of  current  and  emerging  radiation  delivery  systems.  The  phantom  is  easily  adapted  to  use 
in  experiments  with  any  clinical  linac,  tomotherapy  system,  or  even  stereotactic  radiotherapy 
devices.  Ultimately,  as  the  diversity  and  complexity  in  the  technology  of  radiation  therapy  and 
radiobiological  research  increases,  the  biological  and  technological  flexibility  of  the  phantom, 
along  with  its  extreme  efficiency  as  an  experimental  tool,  makes  it  an  immensely  useful  and 
potentially  important  tool  in  the  understanding  of  the  underlying  radiobiology  of  modern 
radiation  therapy. 


5.  Conclusion 

We  have  created  a  simple,  flexible  and  efficient  phantom  for  use  in  three-dimensional  IMRT- 
based  in  vitro  experiments.  The  simultaneous  use  of  multiple  tissue-culture  plates  provides 
a  three-dimensional  grid  of  biological  sample  points,  allowing  for  biological  experiments 
to  be  performed  in  a  realistic  IMRT  setting.  This  provides  a  number  of  improvements  in 
experimental  technique  over  other  previously  used  apparatuses,  including  the  ability  to  sample 
a  large  number  of  points  within  the  dose  delivery  volume  in  a  single  irradiation  session  and  the 
use  of  cell  culture  plates  to  minimize  the  experimental  error.  The  TED  and  film  analysis  show 
that  IMRT-based  dose  distributions  can  be  accurately  and  reliably  projected  into  the  phantom 
under  a  wide  range  of  experimental  setups.  These  factors  make  the  phantom  a  potentially 
important  tool  in  the  understanding  of  radiobiology  principles  that  have  a  direct  impact  on  the 
efficacy  of  radiation  therapy. 
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